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We study the conductivity of granular superconductors in the weak coupling insulating regime. We 
show that it is governed by the hopping of either electrons or Cooper pairs depending on the relation 
between the superconducting gap and the charging energy of a single granule. Local superconducting 
pairing plays an important role in both cases. In particular, in the case of the transport via electron 
hopping the superconducting gap suppresses the inelastic cotunneling processes. We determine 
transport characteristics of an array in different regimes and construct the transport phase diagram. 

PACS numbers: 74.81.-g,74.78.-w,74.45.+c,74.20.-z 



Recent experiments posed the fundamental questions 
concerning mechanisms of conductivity in the insulating 
phase of superconductor granular arrays [l| . The origin 
of the Mott variable range hopping - like resistivity tem- 
perature behavior and the nature of strong enhancement 
of conductivity by the applied magnetic field in the weak 
coupling regime are far from being understood. 

At the same time, significant progress has been re- 
cently achieved in understanding the electronic transport 
in the insulating phase of metallic granular arrays. A long 
standing puzzle of the low temperature resistivity p{T) 
showing stretched exponential temperature dependence, 
Inp ~ T -1 / 2 0, in the insulating phase was explained 
in 0,0,0 as the Mott-Efros-Shklovskii (ES) variable 
range hopping. The key ingredient of the model of Refs. 
[HyLIH is the effect of the electrostatic disorder that lifts 
the Coulomb blockade on a part of sites of a granular 
array providing the necessary low energy electron and 
hole excitations carrying the current by hopping through 
the virtual states of intermediate grains. The specific 
hopping mechanism depends on the temperature range 
0,13: at temperatures T < T\ w 0.1\AE C <5, with E c and 
6 being the charging energy and the mean energy level 
spacing in a single grain respectively, an electron hops 
via elastic cotunneling mechanism such that its energy is 
conserved on each hop. At T > T\ the inelastic cotun- 
neling mechanism Q where the travelling electron cre- 
ates or absorbs the electron-hole excitations dominates 
the transport. The multiple cotunneling as basic mech- 
anism of hopping was recently confirmed in experiments 
on gold nanoparticle multilayers with controlled struc- 
ture and size of the granules [8| . 

Building on these developments we address in this Let- 
ter the low temperature transport in the insulating phase 
of granular superconductors. We show that in this case 
the conductivity is also of the hopping nature but it is 
mediated by hopping of either electrons or Cooper pairs 
depending on the relation between the superconducting 
gap and the charging energy of a single granule. Electro- 
static disorder plays a crucial role in our approach; mod- 
els based on regular superconductor arrays studied exten- 
sively earlier [we refer to the book || for review] show 
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FIG. 1: The diagram shows two regions in the EH ( E c > A ) 
regime where the hopping conductivity goes via elastic (gray) 
and inelastic (white) cotunneling mechanisms in magnetic 
field H vs. temperature T plane. The line H C (T) represents 
dependence of the single grain critical field on temperature. 

activation transport behavior in the insulating phase. 

The immediate effect of superconductivity on the prop- 
erties of granules is two fold: first, the single particle gap, 
A, forms within the each grain. Secondly, the so-called 
parity effect |lfj Hlj | appears: grains that have an even 
number of electrons, N, have the lower energy than those 
that have an "unpaired" odd electron that carries extra 
energy A |"Lq|. We demonstrate that these features af- 
fect strongly transport properties of granular supercon- 
ductors when electrostatic disorder is taken into account. 
Depending on the relation between the local charging en- 
ergy, E c , and the local superconducting gap, A, the low 
temperature transport goes either via hopping of elec- 
trons for E c > A (we hereafter will be referring to this 
mechanism as to EH regime) or via hopping of the Cooper 
pairs for E c < A (CPH regime). We show that this clas- 
sification holds as long as the dispersion in grain sizes can 
be neglected; neither it is affected by the presence of the 
of-diagonal part of the Coulomb interaction. Derivation 
of the conductivity temperature behavior in a general 
case of an arbitrary relation between E c and A is, how- 
ever, cumbersome and below we present the results only 
in the limiting cases corresponding to E c 3> A in EH- 
and to E c <C A in CPH regimes. 

In the EH regime the conductivity, is determined either 
by the elastic or inelastic cotunneling mechanism depend- 
ing on the temperature regime. In the low temperature 
elastic regime, the conductivity is almost unaffected by 
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the gap at E c ^> A and follows the ES law 

a ~ expf-CTo/T) 1 /^ T = &e 2 /£ a ^, (1) 

where e is the electron charge, k is the effective dielectric 
constant of a granular array yjj, is the localization 
length measured in the units of the grain diameter a (we 
consider a model where all grains are of the same size) 
and b is the numerical factor of the order of unity depend- 
ing on the array morphology. The localization length in 
the case of elastic cotunneling is Q 
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where g is of the order of the average tunneling conduc- 
tance; more precisely it is the geometrical average of the 
tunneling conductances g^i+i along the typical tunnel- 
ing path lng= (In g) = (1/JV) J2i=i The energy 
E ~ E c is analogously defined as a geometrical average 
In E — (In E) , where E is the combination of single grain 
electron E\ and hole E % _ excitation energies at A = 



& = 2/{l/E\ + l/Ei). 



(3) 



In the absence of the superconducting gap the elastic 
transport regime crosses over at T w Tj to the regime 
dominated by inelastic cotunneling processes where the 
conductivity is given by ES law with weakly temperature 
dependent localization length 0, 



t£? = 2/ln(E 2 /16nT*g). 



(4) 



In the presence of the gap and at T < A the inelastic 
cotunneling processes are suppressed. Thus the result 
holds only in the interval T > max{A,Ti}. Considering 
the case T c > Ti, we find the new transport regime at 
T < T$ sa ^fj 1 A, which is dominated by inelastic pro- 
cesses but has the activation form: 
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where the typical tunneling order N = \Jbe 2 /16okA ~ 
s/E c / A. Under temperature decease dependence (JjJJ re- 
mains applicable till it matches the elastic result ill 121) at 
temperature T 2 « ^ H A. 

Changing the gap A by varying magnetic field at fixed 
temperature one can drive the system between the elastic 
and inelastic ES regimes traversing the interval T-z(H) < 
T < Ts(H), where conductivity follows the activation 
formula of Eq.© (see Fig.l). This manifests itself as 
a giant negative magnetoresistance in the temperature 
interval (Ti, T3) that exists provided T\ <C T c . Such an 
effect in the weak coupling regime was indeed observed 
in Ref.Q. 

In the CPH regime the conductivity follows the ES law 
with Tq and the localization length given by 




2E--2A 



A-E € 



1/ 7 V v 

i^'i 1 i v 1 N 1 • r 1 1 1 1 • 

-3E C -2E C -E c E c 2E C 3E C -2E C -E c E c 2E € 
a) b) 

FIG. 2: Single particle (thin lines) and two particle (thick 
lines) excitation energies as functions of the applied potential 
in a single grain for the cases a) E c > A and b) E c < A. 
Solid lines represent creation processes (addition of an extra 
electron or a pair) while dashed line represent annihilation 
processes. In both cases the dependencies are 4_B C periodic. 

We see that the localization length decreases with de- 
crease of A, which corresponds to the strong positive 
magnetoresistance. Note further that at gA ~ E c the 
sample transforms into a superconducting statejl^. 

We begin our quantitative analysis with the discussion 
of the ground state charge configuration at T — > 0. First 
we consider an isolated superconducting grain in the pres- 
ence of the random potential V modelling electrostatic 
disorder. The energy of such a grain is 



E = n 2 E c -Vn + P(n + p) A, 



(7) 



where P is the parity function defined as P(2k + 1) = 
1, P(2k) = for an integer k, and n is the number of 
excessive electrons, counted with respect to the number 
of electrons, N°, of a neutral state at V — 0, A = 0. The 
discrete random variable p = for even parity and p = 1 
for odd parity of the charge neutral state N°. 

The effect of the applied potential V on the ground 
state is qualitatively different for different relations be- 
tween E c and A [Tfl HJ. If E c > A (EH), the depen- 
dence of the ground state charge on V has a form of the 
Coulomb staircase with the charge jumps n — > n + I at 
V n = (2n + 1)E C + A cos7r(n + p). In the case E c < A 
(CPH), the dependence of the occupation number on V 
also has a staircase form but electrons jump in pairs, 
n — > n + 2, at the equidistant values V n — E c (2n + 2), 
with even-, n — 2k, for p = - and odd, n — 2k + 1, for 
p = 1, occupation numbers (k is an integer). 

Further, one can see that in the EH regime the Cooper 
pair excitations have a gap, while single particle excita- 
tions are gapless for certain values of V. On the con- 
trary, in the CPH regime electron and hole excitations 
are gapped for all potential values while the pair excita- 
tions are gapless at certain values of V. Shown in Fig.l 
are electron £ + and hole £_ excitation energies 
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£± = (±2n+ l)E c =F V + Acos7r(n +p), 



(8) 
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and pair creation £2+ and annihilation £2- energies 
£ 2 ± =4(±n+l)E c T2V, 



(9) 



as functions of V; the minimal pair excitation energy in 
the EH regime is 2(E C — A) > 0, and the minimal energy 
of an electron excitation in the CPH regime is A— E c > 0. 

We assume that electrostatic disorder is strong taking 
characteristic dispersion of V as Vo > E c . Averaging over 
the potential V results in the finite density of low energy 
electron-hole excitations in the EH regime while in the 
CPH regime only the pair gapless excitations appear. 

The excitation spectrum we have discussed refereed to 
nonintcracting granules. Now we introduce interactions: 



H c = ''</•,'•'", - Vim + P{rn + Pi) A, 



(10) 



where the matrix Efp represents the Coulomb interac- 
tion related to the capacitance matrix Cy as E l J — 
(e 2 /2)C~ j 1 . Remarkably the main conclusion concerning 
the nature of low energy excitations depending on the 
ratio of E c = E" and A obtained for a single grain holds 
for model 1)11) fl . Indeed, the energies of the single and two 
particle excitations on the site i following from (|10fl are 



£1 = 

ci _ 
C 2± ~ 



±fa ~ 

±2pa 



E c + Acos7r(ni 



Pi), 



(11) 
(12) 



where the local potential fa is fa = 2 Y^j E l J n j ~ V%- The 
stability of the ground state charge configuration assumes 
that £±,£2± > that limits possible values of fa. Now 
suppose that fa is such that an electron excitation on 
the grain i is gapless: £\ = 0. The ground state stability 
requires that £ l _, £ 2 +, £\- > giving immediately E c > 
A and £|± > 2(E C — A). This means that gapless electron 
excitations exist only in the case E c > A and that at the 
same time pair excitations have a minimal gap 2(E C — A). 
Analysis of the case £1=0 yields the same result. 

Similarly, for a gapless pair excitation to appear at a 
certain grain i, it is required that £\+ — and stability 
assumes £i, £\, £\_ > 0, meaning that fa + pi has to 
be even and that A > E c , £± > A — E c . This in tern 
assumes that the gapless pair excitations exist only in 
the case E c < A and that single particle excitations has 
the gap A - E c > 0. 

To find the density of low energy excitations we follow 
the standard ES approach For the EH regime we 
require that the energy of the excitation consisting of 
replacing an electron from the site i to the site j were 



positive: £_ 



2E % i > 0, resulting in the 



standard ES expression for the low energy single particle 
density of states (DOS) 

v x (e) = a dl {k/e 2 ) d e d -\ (13) 

where a c i is the dimensionless factor |l5|. In the CPH 
regime we find the density of two particles excitations by 




FIG. 3: The diagrams represent the tunneling probabil- 
ity via elastic (a) and inelastic (b) cotunneling processes. 
The crossed circles stay for the tunneling matrix elements 
t'k k' e^'^' e -l< ^' T ' where phase factors appear due to the 
gauge transformation. Wavy lines represent the average of the 
phase factors (e** (n) e^ (T2) ) with respect to the Coulomb 
action. 

requiring the stability of the ground state with respect to 
a replacement of a pair from site i to the site j: £2-2+ = 
£\_ ~ &E'c > 0. We again arrive at the ES stability 



condition leading to the pair DOS as 

v 2 {s) = a d2 {k/(2e) 2 ) d e d -\ 



(14) 



The two regimes differ only in the doubled charge of the 
Cooper pair, and thus we expect that a d \ ~ 0:^2 • 

Now we turn to the derivation of hopping probabilities. 
We work in the basis of exact eigenstates of noninteract- 
ing isolated grains. Electron tunneling processes between 
the states ki of the grain i and kj of the grain j are rep- 
resented by the elements of the tunneling matrix t l £. k . 
The Coulomb interaction within each grain is accounted 
for via the gauge transformation of the electron fields 
c fc ( T ) Ck (t) e^ r ) , and the phase field 4> appears as a 
renormalization of the hopping matrix elements |l6j 

iy->t lJ e* (T) "* W . (15) 
The phase field cf> is governed by the Coulomb action 

s = - 2? E / dT & + iV *) c « + lV i)- ( 16 ) 

ij 

Without the Coulomb interaction a grain is described by 
the standard Bardeen-Cooper-Schrieffer model. 

The tunneling probability via the elastic process given 
by the diagram a) in Fig. can be written as a prod- 
uct P ~ Y[i 9i,i+i Pi with Pi representing an elementary 
contribution from the grain i 

Pl = iJ d ^ dT ^ G d^)Gd^)e- E ^- E -^. (17) 

Here E\ and E % _ are given by Ea. ((TT|l with A = and 
G{(t) is the superconducting Green function 

1 



G ? (r)=T- f(TEt/T)(l + t/E ( ) 



+f(±Es/T)(l~Z/E e )e E t T \, (18) 

where the upper(lower) sign stays for r > (r < 0), 
f(x) = 1/(1 + e x ) and E± = yj£ 2 + A 2 . Taking the inte- 
grals in (|17fl in the limit T -C A we obtain 



P=S/tt£\ £ l = E l +ttA/2, 



(19) 
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that leads to the localization length of the dependence 
JTJ. Neglecting the term irA/4 in the expression for £ l 
we arrive at the expression J2J for the localization length. 
We would like to note that this correction though being 
small in the considered limit (E c ^> A) can still lead to 
a noticeable negative contribution to the magnetoresis- 
tance in the regime of elastic cotunncling. 

The tunneling probability via inelastic cotunncling 
processes is given by the diagram b) in Fig. J2J that in 
the limit T, A < E c gives 

/oo 
dt[G(-it)G(it)] N e- ist , (20) 
-oo 

where N is the tunneling order, G(t) = f d£ G^(r) and 
e is the electron energy change. The expression (|20f) 
has to be minimized with respect to N under constraint 
NasR/e 2 ~ 1, that follows from the Mott argument for 
finding the typical minimal distance to a site available for 
electron placement within the energy shell e and from ES 
expression for DOS. In the limit N ^> 1 the integral 
in (|20|) can be taken within the saddle point approxima- 
tion leading to ES law 111 II) for temperature T > T3 and 
to the activation behavior J^J at T < T3. 

Deriving the result JHJ) we neglected the possibi lity of 
inelastic cotunneling through the unpaired states [lfj of 
" odd" grains that constitute about a half of all grains in 
the case E c ^> A. One can show that such processes do 
not affect the results since the elementary probability of 
the inelastic cotunncling process via an unpaired state is 
smaller that for the elastic process by the factor of T/E c . 

The tunneling in CPH regime can be described within 
the effective model acting on Cooper pairs 

H = - 2 + (1/2) Y J y ' 

ij i <ij> 

where ip and n — —id/dtp are the Cooper pair phase and 
number operator respectively and J^- = g^-K A/2 |l7| is 
the Josephson coupling between the neighboring grains i 
and j. Tunneling amplitude, A, can be calculated within 
the straightforward perturbation theory in J: 

N 

A~n j w/^ ^2=2/[l/4+ + l/4-], (21) 

where pair excitation energies £\+^2- arc defined by 
Eq. (|P2|) . The tunneling probability, P = A* A, decays 
with the tunneling order N as e~ 2N ^ with the local- 
ization length ^ CPH = 1/ 111(2^2/71- <? A), where E 2 is the 
geometrical average of effective pair excitation energies 
£2 along a typical tunneling path. Since the density of 
states in EH and CPH regimes differ only by the effective 
charge, E% ps 4E leading to Eq.©. 

In conclusion, we have described the hopping trans- 
port in granular superconductors and found that if the 
single grain charging energy exceeds the superconducting 



gap, E c > A, the transport goes via hopping of electrons, 
while in the opposite case, E c < A, hopping of Cooper 
pairs dominates the transport. In the former case we 
predict the negative magnetoresistance, while the latter 
regime exhibits the positive magnetoresistance. We re- 
late the giant negative magnetoresistance observed in the 
insulating phase of the granular superconductors in to 
the suppression of the inelastic cotunncling in the elec- 
tron hopping dominated regime. Transport via Cooper 
pair hopping can be observed in samples with low enough 
grain charging energy, i.e. large enough grains. Such 
regime can also appear as a result of the renormalization 
of the Coulomb energy due to intergranular coupling in 
samples with intermediate coupling strength g ~ 1, in 
particular close to the insulator to superconductor tran- 
sition one expects CP dominated transport. The regime 
g ~ 1 is a subject of our future investigation. 
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